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1.  Introduction 


The  combat  requirements  for  modern  fighter  aircraft  necessitate  that  tiiey  be  highly  maneuverable  and 
able  to  operate  at  high  angle  of  attack.  Since  most  fighter  aircraft  incorporate  a  delta-wing  lifting  surface, 
there  has  been  great  interest  in  studying  the  aerodynamics  of  delta  wings.  'Fhese  flows  are  al.so  rich  in 
fundamental  fluid  mechanics  phenomena,  including  longitudinal  vortices,  vortex  breakdown,  and  shi-ar- 
layer  instability.  Extensive  experimental  and  computational  studies  of  delta  wing  flows  have  been  reported 
in  the  literature.  An  extensive  review  of  computational  techniques  for  modeling  delta-wing  flows  has  been 
given  by  Hoeijrnakers  [1], 

More  recently  there  has  been  an  increasing  interest  in  unsteady  delta-wing  maneuvers,  including  pitching 
to  large  incidence  angle  and  high  rate  rolling  motion.  This  unsteady  flow  environment  introduces  many 
interesting  nonlinear  aerodynamic  phenomena.  The  inherent  unsteadiness  of  the  flowfields  can  also  lead  to 
aerodynamic  behavior  that  may  be  advantageous  or  detrimental  to  aircraft  performance.  Understanding  of 
this  unsteady  aerodynamic  behavior  will  be  vital  in  improving  the  design  and  performance  of  present  and 
future  fighter  aircraft.  Ericsson  and  Reding  [2]  have  provided  an  extensive  review  of  the  fluid  dynamics  of 
tliese  unsteady  separated  flows. 

The  predominant  interest  in  roll  maneuvers  has  been  understanding  the  self-sustained  wing  rook  (die- 
nomena  observed  for  high  sweep  della  wings  at  moderate  to  large  angles  of  attack.  Several  differi'ni 
experiments  in  both  water  and  wind  tunnels  have  investigated  delta-wing  rock  through  flow  visualization 
and  force  and  moment  measurements  [3-7],  Recently,  Arena  [8]  has  also  reported  unsteady  surface  pressure 
measurements  for  a  rocking  delta  wing.  These  experiments  have  provided  valuable  insight  into  the  mech¬ 
anisms  that  lead  to  seif-sustained  wing  rock.  They  are  limited,  however,  by  the  difficulties  encountered  in 
taking  measurements  in  a  dynamic  environment. 

Numerical  simulations  of  delta-wing  rock  provide  an  alternate  means  of  investigating  these  rolling 
motions  and  compliment  experimental  measurements.  Numerical  simulations  allow  for  the  investigation  of 
flow  features  that  may  be  difficult  or  impossible  to  obtain  experimentally.  The  majority  of  computations  for 
rolling  delta  wings  have  assumed  either  conical  (9-11)  or  inviscid  [12]flow  and  their  application  is,  therefore, 
limited,  t’haderjiaii  [13]  has  presented  the  first  computations  for  wing  rock  using  the  three-dimensional, 
thin-layer  Navier-Stokes  equations. 

The  focus  of  the  present  work  is  the  numerical  simulation  of  a  constant  roll  rate,  $  =  0.1325,  maneuver 
from  0"  to  45®.  The  study  of  this  “roll-and-hold”  maneuver  is  useful  in  understanding  the  aerodynamics 
of  a  maneuvering  fighter  aircraft.  While  this  transient  maneuver  does  not  contain  the  full  wing-rork 


I 


motion,  many  of  the  aerodynamic  features  observed  are  also  conifiion  to  wing  rork  1  le  w-  roinjiiit-iJnuis 
therefore,  provide  valuable  information  and  insight  before  proceeding  with  the  more  rosti>  Kiniputa!! 
of  a  complete  rocking  maneuver.  This  work  is  part  of  an  ongoing  study  of  wing  rock  both  with  and  w  ithout 
vortex  breakdown. 
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2.  Governing  Equations 


The  governing  equations  for  the  present  problem  are  the  unsteady,  three-ditnensional,  full  in  ass- aver  age  cl 
Navier-Stokes  equations  written  in  strong  conservation  form  [14]  using  a  general,  tinnvdependent  coordi¬ 
nate  transformation  i;,  t,M: 

f  V  \  dH  _  \  ,dl\.  dG^  dlh.. 

where: 


r  1 

pV 

pV 

pif' 

pu 

^  1 

puU  -t-  C^p 

1 

puV  -f  r;^p 

,,  1 

puW  +CiP 

pv 

•  ^  =  7 

pvU  +CvP 

•  '^  =  7 

pvV  +  i]yp 

’  ^-7 

pvW  -f  CyP 

pw 

pwU  -t-  CjP 

pwV  +  rj^p 

pwW  +  CiP 

pE, 

L  j 

_  pE,U  +  p6^  ^ 

pE,V  +  pV 

pEtW  +  pW 

0 

0 

' 

0 

iij>\ 

^  1 

VxJt\ 

,,  1 

Cr.ni 

Cx.n? 

.  g„  =  7 

lxJt7 

.  H.  =  - 

<rj,2 

Cr,  na 

Vr.Ej 

GJ,3 

IxM 

Cr,6, 

v  =  m  T  »?r«  +  r]yV  +  J],  u:  =  7),  +  V 
=  Cl  +  Cr  n  +  Cy  I'  +  C/  =  Cl  + 


Ef-eG  +  w^) 

A  -  ^ 

0,  _  u,r,j  +  _  l)PrMl  dx.  DG 

i^dui  d^k  ^  2  djt  duj 

dxj  dik  ^  dxi  dik  3^ 

p  =  fii-i-fi, 

xi,X2,X3-^  x,y,z 
Ci.C2.43-^«?.C 
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The  system  of  equations  is  closed  using  the  perfect  gas  law,  Sutherland's  formula  for  viscosifv  and  (hi- 
assumption  of  a  constant  Frandt!  number,  Pr  =  0  72.  Flow  quantities  have  been  tioiidntien.sion.dired  (iv 
their  respective  free  stream  values  except  for  pressure,  which  i.s  nondimen.sionalized  by  twice  the  fre.-stream 
dynamic  pressure,  and  speed  of  sound  which  is  nondimensionalized  by  the  freestream  veiociiy  All  lengths 
have  been  normalized  by  the  root  chord  length  of  the  delta  wing. 
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3.  Numerical  Procedure 


rih>  govcriiinj^  r(ni:iti>>iis  an-  snivi  d  iiuini-ncallv  umii^  tlu'  implicit,  appro.'iimati'ly-fac turcd  algorithm  of 
Beam  and  Warming  [l.'i;  l  lic  equations  are  differenced  using  F.uler  implicit  timc-differencirig  and  second- 
order  accurate  central  differences  for  all  spatial  derivatives  A  blend  of  second  and  fourth  order  nonlinear 
dissipation.  a.s  suggesied  1\\  Jameson  et  al  [Itj],  is  used  to  stabilize  the  central  difference  sclieine  The 
Current  Work  in  whicti  siibs<.)nic  flows  are  investigated,  requires  only  fourth-order  dissipation  In  ilie  |ir<'seiii 
uunierica!  scheme  both  a  full  block  t n-diagcjnal  inversion  schente  and  a  diagonalized  inversion  sclieine  bas<.*d 
<111  till-  diagonal  form  of  tin-  Hearn- Warming  algorithm  developed  by  Pulliam  and  Chaussee  [17]  are  available 

■\  subiti-ration  proceilure  that  luts  been  successfully  used  by  wveral  authors  [18-20]  has  also  been 
iiK'oriiorated  as  an  option  in  the  current  solution  procedure  [21]  The  snbiteralion  strategy  providi’S 
distinct  impruveiuenis  to  the  numerical  scheme  The  stability  limits  of  the  three-factored  algorithm  are 
reia.aed  hy  reducing  the  factorization  error  through  subiterations.  This  greatly  improves  the  efficiency  of 
the  algorithm  Subiteratioiis  also  provid"  a  consistent  means  of  implementing  the  diagonal  form  of  tlie 
ifeam-Waniuiig  algorithin  for  unsteady  flows 

A  fully  \eciori/,ed,  time-accurate  code  luas  be'en  developed  to  implement  the  aforementioned  scheme  [22j. 
fill'  l  oiupiii al  loiial  ri-<|mr>'menis  are  dd  words/gnd- point  anil  a  processing  rate  of  1.8  x  10"*  CPf!  ser/grid- 
poml/iterat ion  fur  the  diagonal  solver  on  a  Cray  !l  I'lie  block  tri-diagonal  solver  requires  appro.viiuaiely 
twin-  the  diagonal  vjlver  processing  rate  [  fie  code  lias  been  validated  for  a  variety  of  steady  and  unsteady 
flows  Both  su|)ersoiiic  and  subsonic  flow,s  over  della  wing.s,  including  vortex  breakdown  (21,23-25)  have 
l)|■^•n  !  <uupuieil  with  this  code  I  iisieady  flow  simulations  have  been  carried  out  for  both  a  pitching  slender 
body  of  revolution  [2hj  and  a  (utrhing  delta  wing  [25).  Finally,  the  code  has  been  usea  to  calculate  steady 
and  luisleady  licjrsi'shoe  vortex  flows  occurring  at  a  cyiinder/flat  plate  juncture  [22,27],  These  results  have 
pri/vnh-d  a  liroad  vainlation  bav  for  both  steady  and  unsteady  flows. 

For  tie-  i>ody  motion  to  be  considered  in  tins  study  fi  e.,  delta- wing  roll),  a  general  time-dependent 
<  oordiiiat<  irausforiiiation  is  introduced  The  rigid,  nondeforming  grid  is  allowed  to  roll  with  the  body  and 
tin-  grill  mi.tioii  IS  treated  throiigli  llie  general  coordinate  transformation.  This  approach  eliminates  tin- 
III  I  *!  tor  gi-ii'  rating  mulli|>le  grids  Also,  siiicr?  the  equations  are  written  in  an  inertial  frame  of  reference, 
no  .idditu.iial  terms  are  r*-<piir<’d  in  tin-  governing  e<|uations. 


4.  Grid  Structure  and  Boundary  Conditions 


The  delta-wing  configuration  used  for  tlie  present  compulations  corresponds  to  llie  delia-winj’  model 
used  by  Arena  [8j  in  his  experimental  wing  rock  studies.  The  wing  has  an  80°  leading  edge  swee,>,  a 
root  chord  of  16  5/8  in,  a  thickness  of  1/4  in  and  a  45°  bevel  at  the  leading  edge.  In  the  computational 
configuration  the  trailing  edge  has  also  been  beveled  (9.3°)  to  provide  for  a  smooth  closure  of  the  wing. 

The  grid  chosen  for  the  current  study  h^  an  H-0  structure,  with  O-grids  being  stacked  axially  along 
the  delta  wing  (Fig.  1).  This  grid  topology  provides  good  resolution  of  the  leading  edge  vortices  and  the 
apex  region  of  the  delta  wing.  The  baseline  grid  used  for  these  computations  consists  of  219  points  around 
the  body,  81  points  normal  to  the  body  with  a  minimum  spacing  at  the  wall  of  il;  =  O.OOOl  and  89  points 
in  the  axial  direction  with  a  constant  grid  spacing  of  A-r  =  0.025  over  the  predominant  portion  of  the  wing 
The  grid  extends  1  chord  upstream  of  the  wing,  1.5  chords  downstream  of  the  wing  and  1  ■')  chords  in  the 
normal  direction  from  the  wing. 

The  boundary  conditions  for  the  delta-wing  roll  problem  are  implemented  as  follows  Characteristic 
boundary  conditions  [28]  are  applied  at  the  far  field  and  upstream  boundaries.  At  the  downstrf  am  bound¬ 
ary  first  order  extrapolation  of  the  interior  flow  values  is  used.  On  the  delta-wing  surface  the  following 
conditions  are  applied: 

=  Ch 

=  0 

=  •  n 


U 

dri 

dn 


where  Vf,  and  aj  are  the  velocity  and  acceleration  of  the  point  on  the  body  defined  as: 

Ui  =  J  X  (fj  -  r„) 

dQ 

04  =  -37  X  (fi  -  To)  -1-  w  X  [w  X  (fj,  -  r„)] 

at 

and  n  is  the  surface  normal  vector  (see  Fig.  2).  For  the  case  of  roll  about  a  fixed  axis  considered  here, 
r„  =  0  and  w  =  <t>i  where  4>  is  the  roli  rate. 

For  the  H-0  grid  structure  adopted  for  the  present  problem  a  singular  line  emanates  from  the  apex  of 
the  delta  wing  and  a  singular  line  and  plane  are  located  downstream  of  the  trailing  edge.  These  singular 
surfaces  are  the  continuation  of  the  grid  from  the  wing  surface  ahead  and  downstream  of  the  delta  wing. 
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The  flow  variables  on  tliese  surfaces  are  obtained  by  extrapolating  the  conserved  variables  to  the  surface 
and  tlien  performing  t!)e  appropriate  averaging  of  these  values. 


9 


5.  Results 


5.1  Fixed  Roll  Angle 

Two  fixed  roll  angle  cases  =  0“  and  ^  =  45®)  are  considered  to  deiiionslrate  tlio  capaljilitn’s  of 
the  numerical  scheme  before  computing  the  dynamic  roll  maneuver.  The  flow  conditions  are  a  Reynolds 
number  of  400,000,  freestream  Mach  number  of  0.2  and  an  angle  of  attack  of  30®.  These  conditions 
correspond  to  the  experiments  of  Arena  [8].  These  experimental  measurements  are  used  as  a  validation 
base  for  comparison  with  the  present  computations  for  rolled  delta  wings. 

Computations  for  both  laminar  and  turbulent  flow  cases  were  performed.  For  the  turbulent  flow  cases 
the  Baldwin-Lomax  turbulence  model  [29]  is  used.  The  turbulence  model  has  been  modified  similarly  to 
Ref.  [30]  to  account  for  the  large  vortical  structures  over  the  delta  wing  by  restricting  the  search  for  the 
peak  in  the  Baldwin-Lomax  F-function  to  a  normal  distance  from  the  body  of  no  more  than  30%  of  the 
local  wing  span.  The  flow  is  assumed  to  be  turbulent  over  the  whole  extent  of  the  wing 

I’he  st'usitivity  of  the  solution  to  grid  resolution  and  boundary  location  is  also  investigated.  Compu¬ 
tations  are  performed  on  a  second  grid  containing  299  x  81  x  89  grid  points.  Both  the  baseline  grid  and 
the  refined  grid  have  the  same  axial  distribution  of  grid  points  with  the  baseline  grid  having  100  points 
and  the  refined  grid  140  points  distributed  spanwise  across  the  upper  surface.  The  grid  distribution  in  the 
body  normal  direction  on  the  refined  grid  has  been  modified  to  enhance  the  resolution  in  the  region  of  the 
vortex  cores. 

Figure  3  compares  the  computed  surface  pressure  at  several  axial  locations  with  the  experimental 
measurements  of  Arena  [8]  for  a  0®  roll  angle.  The  computed  solution  agrees  well  with  the  experiments, 
capturing  all  significant  features  observed.  At  X/L  =  0  9  the  minimum  pressure  peak  is  slightly  under- 
predicted.  The  discrepancy  in  the  lower  surface  pressure  at  X/L  =  0.9  is  due  in  part  to  the  negative 
camber  introduced  by  the  lower  surface  trailing-edge  bevel  used  in  the  computational  model.  Any  other 
discrepancy  between  the  experimental  and  computed  pressures  is  well  within  the  experimental  uncertainty 
XCp  0.079—  0.084  [8],  The  effect  of  grid  boundary  location  on  the  computed  solution  has  also  been 
investigated  for  this  case.  The  distance  from  the  wing  to  all  boundaries  was  doubled  with  no  noticeable 
changes  observed  in  the  solution. 

The  next  case  considered  is  a  delta  wing  at  a  fixed  roll  angle  of  45®.  The  vortex  structure  over  the  rolled 
delta  wing  is  seen  in  Fig.  4.  The  right  vortex  (downward  leading  edge)  has  moved  inboard  and  towards  Ihi' 
surface  of  the  wing  creating  an  extensive  secondary  flow  region.  The  left  vortex  (upward  leading  edge)  has 
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Figure  4;  (Contours  of  the  Axial  (.Component  of  Vorticity;A7i  =  0.6.0  =  45® 


moved  outboard  of  the  wing  and  away  from  the  surface  and  its  iidlucnce  on  the  flow  over  the  upper  surface 
of  the  delta  wing  is  greatly  diminished.  The  shear  layer  that  rolls  up  to  form  this  vortex  now  separates 
from  the  lower  surface  bevel  rather  than  the  upward  leading  edge. 

Figure  5  shows  a  comparison  of  the  computed  surface  pressure  at  XjL  -  0.6  for  both  the  base  grid 
and  the  refined  grid  with  experiment  [8].  Both  computed  solutions  agree  well  with  the  experimental 
measurements,  with  only  a  slight  overexpansion  in  the  region  under  the  right  vortex.  Improvements  in  the 
solution  in  this  region  are  seen  for  the  refined  grid.  Elsewhere  the  two  solutions  are  virtually  identical.  This 
figure  also  reflects  the  reduction  in  influence  of  the  left  vortex  on  the  surface  pressure  with  the  expansion 
peak  due  to  the  left  vortex  being  eliminated. 

Figure  6  compares  the  experimental  surface  oil  flow  with  the  simulated  oil  flow  patterns  for  the  com¬ 
putations  on  the  baseline  and  refined  meslu's.  flood  (lualitative  agreement  between  the  main  fi'atures  of 
tJie  experimental  and  simulated  oil  flow  patterns  is  obtained.  Both  the  secondary  and  tertiary  separation 
lines  and  the  attachment  of  the  downward  leading  edge  vortex  on  the  upper  surface  are  captured.  The 
curvature  of  the  separation  lines  as  they  approach  the  trailing  edge  observed  in  the  experiment  is  also  seen 
in  the  computed  oil  flow  patterns.  Some  differences  in  the  aft  portion  of  the  wing  between  the  tertiary 
separation  line  and  the  leading  edge  are  seen,  but  the  true  flow  topology  in  this  region  is  not  discernible 
even  in  the  experimental  oil  flow  picture. 

Since  there  is  not  a  unique  relationship  between  the  surface  topology  and  the  flow  pattern  above  the 
wing,  the  cross-flow  topology  is  shown  in  Fig.  7  at  an  axial  location  XfL  -  0.6.  In  this  figure  features 
are  marked  as  follows:  S  are  saddle  points,  F  are  foci  and  N  are  nodes.  Primed  quantities  denote  half- 
saddles  or  nodes.  The  subscripts  A  and  S  denote  attachment  and  separation,  respectively.  The  +  and 
-  superscripts  denote  stable  (inward  spiraling)  or  unstable  (outward  spiraling)  foci,  respectively.  This 
cross-flow  topology  sati.sfies  the  topological  rule  for  a  two-dimensional  plane  cutting  a  tlutH'-dimeiisional 
body  [.‘11]  (53^+|5Z/v')  ~  (IZ.S' +  j  IT.?' )  =  -1.  The  flow  topology  at  this  location  on  the  delta  wing 
differs,  however,  from  the  standard  topology  commonly  described  in  the  lileratnre.  No  focus  i.s  .•v.s.socialed 
with  the  left  vortex  which  has  moved  outboard  and  away  from  the  wing.  An  unstable  focus  with  a  stable 
limit  cycle  is  associated  with  the  right  vortex.  This  cross-flow  topology  has  also  been  observed  for  the 
0°  roll  angle  ceise.  Furthermore,  Visbal  [25]  has  found  limit  cycle  behavior  in  the  cross-flow  topology  for 
pitching  delta  wings.  He  provides  a  more  detailed  discussion  of  cross-flow  topology  of  delta  wings  including 
the  occurrence  of  limit  cycles  and  stable  and  unstable  foci.  In  the  secondary  and  tertiary  flow  regions, 
all  the  foci  are  unstable  foci.  The  stable  focus  located  above  the  delta  wing  and  near  the  leading  edge  is 
associated  with  the  accumulation  of  vorticity  in  the  shear  layer  emanating  from  the  leading  edge  of  tlie 
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delta  wing  (see  Fig.  4). 

Turbulent  flow  solutions  have  also  been  computed  for  roll  angles  of  0®  and  45®.  Figure  8  compares  the 
Computed  surface  pressure  coefficients  for  turbulent  flow  with  the  laminar  flow  results  and  the  experimental 
measurements  for  the  45®  roil  case.  4'lie  turbulent  flow  solution  shows  a  higher  suction  under  the  primary 
Vortex  and  a  reduction  in  the  influence  of  the  secondary  flow.  Tliis  behavior  of  the  surface  pressure  i.-; 
consistent  with  the  experimental  ntea.suremenl.s  of  Hummel  [32]  for  a  76®  sweej)  delta  wing  at  0®  roll  aitgle 
where  the  boundary  layer  was  artificially  tripped  halfway  between  the  leading  edge  and  the  centerline  of 
the  wing.  The  comparison  with  the  experiments  of  Arena  [8j  also  further  confirms  that  the  experimental 
flow  is  laminar. 

Table  5.1  gives  a  comparison  of  the  computed  vortex  core  locations  v’.th  the  experimental  core  locations 
at  XI L  =  0.6.  The  experimental  measurements  are  subject  to  an  uncertainty  of  ±0.083.  The  vortex  core 
locations  for  all  cases  lie  well  within  the  experimental  error  range.  The  turbulent  vortex  cor-?s  are  further 
outboard  of  the  corresponding  laminar  flow  core  locations.  This  is  due  to  the  reduced  size  of  the  secondary 
vortex  for  turbulent  flows. 

5.2  Dynamic  Roll  Maneuver 

Figure  9  shows  the  roll  maneuver  to  be  computed.  The  maneuver  consists  of  a  constant  rate  roll  from  0®  to 
45°  with  a  smooth  acceleration  and  deceleration.  The  nondimensional  roll  rate  is  <1>  =  0  1325.  This  is  the 
mean  roll  rate  for  one  cycle  in  the  wing-rock  experiments  of  Arena  [8].  Since  no  experimental  data  exist 
for  this  maneuver,  the  computations  were  run  for  a  much  lower  Reynolds  number  than  the  static  cases, 
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Re  —  10,000.  This  provides  the  most  confidence  m  the  computed  results  by  elimiimltn)j,  any  q^ll^^l()lls 
of  turbulence,  transition,  or  the  suscofilibility  to  and  effect  of  .shear-layer  instabilities  |2I)  Kiirtln  riiiore, 
computing  at  lower  Reynold.s  numbers  provides  improved  grid  resolution  for  tlie  same  iiuiub.u  i  I  grid 
points. 

Figures  10  and  11  present  the  effects  of  lower  Reynolds  number  on  the  flow  Uit  o  1.'/  roll  angle 
Figure  10  indicates  the  effect  on  the  surface  pressure  at  Xj'L  —  0.6.  Only  a  very  small  effect  is  seen  on  the 
lower  surface.  On  the  upper  surface,  the  lower  Reynolds  number  ca.se  has  significaiilly  less  expansion  nn  lei 
the  right  vortex.  The  peak  minimum  pressure  for  the  lower  Reynolds  number  care  al.so  occurs  inboard  of 
the  peak  for  the  nigher  Reynolds  nurnher  case.  These  difFcrence.s  in  the  surface  pre.ssure  may  be  .attributed 
to  two  causes.  The  right  vortex  is  stronger  for  the  higher  Reynolds  number.  -  -3.9(5.  than  for  the 

lower  Reynolds  number,  —  —1.86.  The  .second  cau.se  i.s  a  difference  in  the  secondary  How  f  at  lire's 

Figure  11  niay  be  compared  with  F'ig  6  to  determine  the  effect  of  lower  Reynolds  iiiimlier  on  the-  surface 
oilflow  pattern.  Distinct  differences  are  ob.served  in  the  secondary  flow  region  with  the  lower  Reynolds 
number  case  showing  no  tertiary  separation.  Fxamination  of  the  cro.s.s-  flow  topology  revealed  that  one 
large  secondary  vortical  structure  exists  for  (he  low  Reynolds  number  case  This  larger  secoiidary  vortex 
caii.ses  llie  right  vortex  to  he  located  further  inboard  thus  shifting  the  peak  in  minimum  pressure.  '|  he  low 
Reynolds  number  case  also  ha.s  a  more  e.xten.sive  rever.se  flow  region  towards  iJie  trailing  edge.  Kven  with 
the  differences  that  have  been  noted,  the  low  Reynolds  number  computations  for  a  forced  rolling  motion 
should  provide  valuable  insight  into  the  overall  dynamics  and  behavior  of  the  primary  leading  edge  vortices. 

The  computation  of  the  roll  maneuver  for  the  delta  wing  is  very  costly,  requiring  approximately  ‘225  CPU 
hours  on  a  Cray  II.  In  order  to  minimize  the  computing  time  required,  the  diagonal  form  of  the  algorithm 
is  used  with  two  subiterations.  The  subiterations  reduce  the  effect  of  tlie  temporal  errors  introduced  in 
the  diagonalization  proccs.s.  This  technique  for  computing  unsteady  flow  fields  due  to  dynamic  motion 
has  been  compared  with  the  use  of  a  fully  time  accurate  block  tri-diagonal  solver  for  the  problem  of  2-D 
dynamic  stall.  A  comparison  (not  included)  up  to  and  through  the  formation  of  the  dynamic  stall  vortex 
.showed  virtually  no  difTereiice  in  tlie  two  solutions.  The  lime  step  used  for  the  unsteady  computation  is 
Af  =  0.001. 

The  thin-layer  Navier-Stokes  equations  are  also  solveil  in  order  to  rediiri-  the  compiitalioii  time.  VVcIi- 
ster  and  Sliang  (25)  have  shown  that  the  thin-layer  and  full  Navier-Stokes  equations  give  nearly  identical 
.solutions  for  delta  wings  at  modest  angles  of  attack.  For  the  case  of  a  dynamic  motion.  Stanek  and  Vis- 
bal  (26]  have  shown  that  no  aiipreciahle  difference  exists  between  the  thin-layer  and  full  Navier-Stokes 
.solutions  for  a  pitching  ogive.  In  the  present  computations  the  vi.scous  effects  are  properly  re.solved  in 
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regions  near  the  wing  aiul  where  the  shear  layer  that  fdVil.s  the  |)riiiiary  vorticd's  is  foriiied.  In  both  these 
regions  the  tliin-layer  approximation  retains  the  relevant  viscous  terms. 

Figure  12  shows  the  typical  behavior  of  the  vortices  in  a  cross-flow  plane  during  the  rolling  motion. 
Contours  of  the  axial  component  of  vorticity  at  X/L  —  0.9  are  shown  for  various  roll  angles  during  the 
maneuver.  Also  included  in  this  figure  is  the  corresponding  surface  pressure  coefficient  at  X/L  =  0.9.  The 
vorticity  contours  show  the  general  dynamic  behavior  of  the  vortices  during  the  roil  maneuver.  As  the 
whig  rolls,  the  right  vortex  (downward  leading  edge)  moves  inboard  on  the  wing  and  toward.s  the  surface. 
The  left  vortex  moves  ontboar<l  of  the  wing  and  away  from  the  surface.  Up  to  approximately  4>  —  ilO® 
the  right  vortex  moves  inboard  and  towards  the  wing  and  the  secondary  flow  region  underneath  the  right 
vortex  grows  in  size.  Subsequently,  the  secondary  flow  region  is  forced  towarils  the  wall  and  diminishes  in 
size.  When  the  wing  stops,  a  new  secondary  flow  region  emerges  with  the  predominant  portion  being  well 
inboard  of  the  leading  edge.  As  the  left  vortex  moves  outboard  and  away  from  the  surface,  the  secondary 
flow  a.ssociated  with  this  vortex  decreases  in  strength  and  eventually  disappears. 

The  actual  motion  of  the  vortex  cores  is  plotted  in  Fig.  13.  The  vortex  core  location  was  defined  as  the 
point  of  minimum  pressure  in  the  core  region.  After  the  initial  acceleration,  the  vortices  move  in  basically 
a  linear  fashion  until  the  wing  stops.  When  the  wing  stops,  the  left  vortex  shows  a  rapid  increase  in  the  Z 
(body  normal)  direction  to  its  final  position.  This  time  lag  in  the  normal  position  of  the  left  vortex  core 
is  similar  to  the  lag  observed  by  Arena  and  Nelson  [33]  for  the  vortex  coming  from  the  upward  moving 
edge  of  a  rocking  delta  wing.  The  change  in  position  of  the  right  vortex  after  the  wing  stops  is  mainly 
attributable  to  the  reemergence  of  the  secondary  flow  region. 

An  indication  of  the  vortex  strength  may  be  obtained  by  examining  the  minimum  pressure  in  the  vortex 
core.  Figure  14  gives  the  variation  of  the  strength  of  the  vortex  cores  as  the  wing  rolls  by  plotting  the 
minimum  pressure  coefficient  in  each  core  at  XfL  —  0.9.  Up  to  <  ss  2.0  (<i>  as  14.6®)  the  strength  of  the 
right  vortex  increases.  After  this  time  the  strength  drops  off  until  t  as  7.45  (1.275  characteristic  times  after 
th<'  wing  stops)  where  the  strength  level.s  off  to  its  final  value.  'I'he  left  vorti'x  consisU'iitly  losi's  strength 
imlil  I  as  7.15.  Similar  behavior  is  observed  if  either  the  maximum  total  pressure  loss  or  axial  (n)  velocity 
in  the  core  is  used  as  indication  of  the  vortex  strength. 

The  behavior  of  the  vortices  during  the  roll  maneuver  may  be  simply  understood  by  considering  the 
effective  angle  of  attack  and  sideslip  angle  (or  alternatively  the  effective  sweep  angle)  during  a  roll  maneuver. 
From  purely  trigonometric  considerations  (no  roll  motion  considered)  the  effective  angle  of  attack  and  the 
effective  sideslip  angle  are  given  by; 
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Figure  12;  Vortex  Dynamics  during  the  Roll  Maneuver:^/ L  =  0  9 
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Figure  l;$:  Core  [.oration  Motion  at  A’/A  =  O.i) 
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F  igure  14:  Vortex  Strength  during  Roll  at  A'/i  =  0.9 
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Figure  15  shows  ««//  /?e//  for  Ih®  30®  angle  of  attack  case  being  considered.  Veerhaageii  and  Naard- 

ing  [34]  have  shown  that  for  delta  wings  the  windward  leading  edge  vortex  moves  inward  and  closer  to 
the  surface  of  the  wing  while  the  leeward  vortex  moves  outward  and  away  from  the  wing  for  increasing 
sideslip.  This  is  consistent  with  the  vortex  motion  observed  for  the  rolling  delta  wing  Vc-erhagen  and 
Naarding  also  show  that  as  the  sideslip  is  increased  the  strength  of  the  windward  vortex  increases  while 
the  stre.  gth  of  the  leeward  vortex  decreases.  Furthermore,  it  is  known  that  the  strength  of  the  vortices  over 
a  delta  wing  decreases  with  decreasing  angle  of  attack.  Using  these  observations  and  Fig.  15,  the  behavior 
of  the  strength  of  the  vortices  during  the  roll  maneuver  may  be  explained.  The  rolling  of  the  delta  wing 
leads  to  a  reduction  in  the  effective  angle  of  attack  and  a  positive  sideslip  angle.  The  net  effect  on  tlie 
left  vortex  (leeward)  i.s  a  continuous  reduction  in  the  strength  of  the  vortex  as  seen  in  Fig.  11.  Figure  15 
shows  that  for  small  roll  angles  (up  to  approximately  20®)  the  effective  angle  of  attack  reduces  very  slowly 
while  the  effective  sideslip  angle  increases  rapidly.  As  higher  roll  angles  are  achieved  the  effective  angle  of 
attack  drops  much  more  rapidly.  This  leads  to  an  initial  increase  in  the  strength  of  the  right  vortex  as 
the  wing  rolls  to  approximately  17®  since  the  effect  of  increasing  sideslip  dominates.  After  a  roll  angle  of 
approximately  22®  the  strength  of  the  right  vortex  starts  to  rapidly  drop  as  the  effect  of  the  reduction  in 
effective  angle  of  attack  increases  and  the  increase  in  the  strength  of  the  vortex  due  to  increasing  sideslip 
levels  off. 

The  effect  on  the  surface  pressure  of  the  variation  in  vortex  position  and  strength  can  he  seen  by 
referring  to  Fig.  12.  The  behavior  of  the  surface  pressure  at  X/L  =  0.9  is  typical  of  other  cros.s-plane 
locations.  As  the  delta  wing  rolls,  the  large  expansion  due  to  the  left  vortex  continually  drops  as  the 
inllueuce  of  this  vortex  on  the  surface  is  reduced  and  it  reduces  in  strength.  The  expansion  peak  due  to  the 
rigid  vortex  initially  moves  inboard  and  grows  until  t  =  2.65  {<j>  =  19.2°)  where  the  value  of  this  peak  then 
starts  to  diminish.  This  correlates  with  the  point  where  the  right  vortex  starts  to  rapidly  lose  strength  due 
to  the  roll  angle.  During  the  roll  maneuver  the  expansion  region  due  to  the  right  vortex  covers  a  growing 
portion  of  the  wing  as  this  vortex  moves  inboard. 

The  effect  of  the  variation  in  surface  pressure  on  the  lift  and  roll  moment  coefficients  is  seen  in  Figs.  16 
and  17.  Figure  16  shows  that  the  lift  coefficient  drops  rapidly  as  the  wing  rolls,  dropping  by  44%  of  its 
original  value  when  the  wing  stops.  The  tilting  of  the  normal  force  vector  as  the  wing  rolls  provides  a 
significant  portion  of  the  loss  in  lift.  The  behavior  of  the  roll  moment  coefficient  during  the  roll  maneuver 
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Figure  15:  Effective  Angle  of  Attack  aiul  Sideslip 
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Figure  17;  Roll  Moment.  Coefficient. 
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is  seen  in  Fig.  17.  The  initial  acceleration  of  the  wing  rreat«'s  a  sharp  rise  iti  (he  roll  niortir  iit  .As  tin  wsag 
rolls  at  a  constant  rate  the  roll  moment  decreases  until  it  reaches  a  niininmm.  Cl  -  -0  0104  at  f  ^  t 
(o  =  The  roll  moment  then  increa.ses  until  it  drops  sharply  when  the  wing  <le(e|praies  to  a  stop 

This  elTect  on  the  roll  moment  coefficient  can  he  explained  hy  considering  tlie  surface  pn  ssuri-,  f  ig  12 
During  the  initial  part  of  the  roll,  the  developiiieiil  of  the  roll  rnomeiil  is  dominated  by  the  decrease  (;f 
the  favorable  roll  moment  due  to  the  left  vortex  and  the  increase  in  the  opposing  monn  til  due  to  th<‘  right 
vortex.  As  the  wing  approaches  and  rolls  through  the  30*  roll  angle,  the  favorable  roll  moment  from  the 
left  vorle.x  has  reached  a  minimum  and  the  opposing  moment  from  the  right  vortex  has  started  to  level  ofT 
The  fight  vortex  has  now  moved  significantly  inboard  and  starts  to  influence  the  pressure  on  the  left-half 
of  the  wing,  enhancing  the  favorable  roll  moment  from  this  half  of  the  wing.  This  combination  leads  to  the 
development  of  the  minimum  in  the  roll  moment  coefficient. 

Both  the  lift  coefficient  and  roll  moment  coefficient  show  long  relaxation  times  to  attain  their  final 
values,  approximately  3.0  characteristic  times.  Both  primary  vortices  have  achieved  their  final  strength 
and  position  hy  1.3  characteristic  times  after  the  wing  motion  has  stopped.  The  additional  relaxation  time 
reipiired  for  the  lift  and  roll  moment  is  mainly  attributable  to  the  time  required  for  ♦he  secondary  flow 
to  establish  itself  after  the  wing  stops.  The  low  Reynolds  number  for  these  computations  enhances  tfii.s 
secondary  (low  effect. 

Finally,  the  dynamic  and  static  data  are  compared  for  a  roll  angle  0  =  26*  to  ascertain  dynamic  effects 
Talile  5,2  compares  (he  vortex  locations,  vortex  strength,  lift  coefficient  and  roll  moment  coefficient,  I'lie 
dynamic  location  of  the  right  vortex  is  outboard  and  closer  to  the  wing  than  the  static  location.  The 
dynamic  location  of  the  left  vortex  is  inboard  and  closer  to  the  wing  than  the  static  location.  (See  also 
Fig.  13.)  The  vortex  strength  is  greater  for  the  dynamic  vortices.  Figure  18  show-s  the  corresponding  effect 
on  the  surface  pre.ssure.  Only  a  small  dynamic  effect  is  .seen  on  the  lower  surface  with  a  more  significant 
effect  apparent  on  ilie  upper  surface  pre.ssures.  The  difference  in  the  upper  surface  pressure  distribution 
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is  consistent  with  the  differences  in  vortex  location  and  vortex  strengtii  previously  noted.  Cioo  .lly.  tin- 
lift  coefficient  shows  little  difference  between  the  dynamic  and  static  coiidiiion,  I'he  static  roll  inornen' 
coefficient,  however,  shows  a  larger  restoring  moment  than  the  dynamic  case.  This  is  primarily  dui-  to  tlx 
lag  in  the  position  and  strength  of  the  left  vortex  during  the  roll  maneuver. 
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6.  Summary 


Computations  of  a  constant  roll  rate,  =  0.1325,  maneuver  from  0®  to  45®  using  a  thrw-iJinii  ii.sional, 
Beam- Warming,  Navier-Stokes  solver  have  been  presented.  To  demonstrate  the  capabilities  of  the  nuinericnl 
procedure,  computational  results  for  fixed  roll  angles  of  0®  and  45®  and  a  Reynolds  number  of  400,000  are 
compared  with  the  experimental  measurements  of  Arena  [8].  Good  agreement  with  the  measured  surface 
pressure  coefficients  is  obtained  for  both  roll  angles.  The  simulated  surface  oil  flow  showed  e.xcelieiii 
qualitative  agreement  with  the  experimental  surface  oil  flow. 

The  computations  of  the  roll  maneuver  were  performed  for  a  significantly  lower  Reynohls  numlier.  Rt  ~ 
10,000,  since  no  experimental  measurements  were  available  for  thus  maneuver.  'I’his  provides  for  added 
confidence  in  the  computed  results  by  eliminating  questions  of  turbiihmce',  transition,  and  the  .susce|>tihility 
to  .shear-layer  instabilities.  The  reduction  in  Reynolds  number  primarily  effects  the  secondary  flow  region 
and  reduces  the  strength  of  the  primary  vortices. 

The  dynamical  behavior  of  the  vortices  at  a  typical  cross  plane  was  described  for  the  roll  maiieuvet. 
The  right  vortex  (downward  leading  edge)  moves  inboard  and  towards  the  surface  while  the  left  vortex 
(upsvard  leading  edge)  moves  outboard  and  away  from  the  surface.  A  lag  in  the  body-normal  position  of 
the  left  vortex  similar  to  the  lag  observed  for  delta  wing  rock  was  noted.  The  left  vortex  continually  loses 
strength  during  the  roll  maneuver.  The  right  vortex  initially  gains  strength  but  then  rapidly  loses  strength 
as  higher  roll  angles  are  achieved.  A  simple,  quasi-static  explanation  of  this  vortex  beha.ior  based  on  the 
effective  angle  of  attack  and  sideslip  angles  during  roll  is  given. 

The  effect  of  the  vortex  behavior  on  the  lift  coefficient,  roll  moment  coefficient,  and  surface  pressure 
during  the  roll  maneuver  has  also  been  described.  The  lift  coefficient  drops  to  44%  of  its  original  value 
at  -  45®.  The  roll  moment  coefficient  develops  a  minimum  at  d*  =  31.55®.  The  developmeni  v>f  thi" 
minimum  can  be  explained  by  the  dynamic  behavior  of  the  vortices  and  their  corresponding  effect  on  the 
surface  pressure  coefficient. 
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